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We report the measurement of muons and muon-induced phosphorescence in DM-Icel7, a Nal(Tl) 
direct detection dark matter experiment at the South Pole. Muon interactions in the crystal are 
identified by their observed pulse shape and large energy depositions. The measured muon rate in 
DM-Icel7 is 2.93 ± 0.04 [j./crystal/day with a modulation amplitude of 12.3 ±1.7%, consistent with 
expectation. Following muon interactions, we observe long-lived phosphorescence in the Nal(Tl) 
crystals with a decay time of 5.5 ± 0.5 s. The prompt energy deposited by a muon is correlated 
to the amount of delayed phosphorescence, the brightest of which consist of tens of millions of 
photons. These photons are distributed over tens of seconds with a rate and arrival timing that do 
not mimic a scintillation signal above 2keVee. While the properties of phosphorescence vary among 
individual crystals, the annually-modulating signal observed by DAMA cannot be accounted for by 
phosphorescence with the characteristics observed in DM-Icel7. 

PACS numbers: 95.35.+d, 29.40.Me, 78.60.-b, 95.85.Ry 

Keywords: DM-lce, sodium iodide, muons, phosphorescence, annual modulation, dark matter, direct detec¬ 
tion, WIMP, South Pole, IceCube 


I. INTRODUCTION 

Astronomical evidence converges on a 27% cold dark 
matter component of the universe [1] , but the properties 
of dark matter remain largely unknown. The Weakly In¬ 
teracting Massive Particle (WIMP) is a theoretically mo¬ 
tivated dark matter candidate [2, 3], and a number of di¬ 
rect detection dark matter experiments are underway [4] 
to search for evidence of WIMP-nucleon scattering [5, 6]. 

The WIMP interaction rate in terrestrial-based detec¬ 
tors is expected to have an annual modulation due to the 
Earth’s motion relative to the galactic WIMP halo [7, 8] . 
The only claim of dark matter detection comes from the 
DAMA/Nal and DAM A/LIBRA experiments at the Lab- 
oratori Nazionali del Gran Sasso (LNGS) in Italy. Com¬ 
bined, these Nal(Tl) detectors have observed a 9.3cr sig¬ 
nal over 14 annual cycles [9], while other experiments 
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have reported null results in conflict with the observed 
modulation [10-21]. Several experiments aim to directly 
test DAMA’s dark matter claim using the same target 
material of Nal(Tl) [22-26]. 

A number of hypotheses for the origin of DAMA’s an¬ 
nually modulating signal have been proposed, ranging 
from alternative dark matter models [27] to potential 
sources of background [28-31]. No background has yet 
been shown to completely account for the DAMA mod¬ 
ulation signal [32-37]. Atmospheric muons are one such 
background, the production of which in the Earth’s at¬ 
mosphere is correlated with the air density in the up¬ 
per atmosphere and therefore exhibits a seasonal modu¬ 
lation [38-43]. In this paper we report the measurement 
of the muon rate and characterization of muon-induced 
phosphorescence in DM-Icel7, a NaI(Tl)-based direct de¬ 
tection dark matter experiment currently operating at 
the South Pole [44]. 


II. DM-ICE17 DETECTOR 

DM-Icel7 consists of two identically-designed detec¬ 
tors, Det-1 and Det-2, deployed 2457m (2200m.w.e.) 
deep in the South Pole ice [44]. The Nal(Tl) crystals 
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are located 7.5 m below the bottom detectors of the Ice- 
Cube Neutrino Observatory (IceCube) [45-49]. Det-1 is 
located near the center of the IceCube array while Det-2 
is located on the array edge, as shown in Figure 1. The 
glacial ice provides a high radiopurity environment, sta¬ 
ble running conditions, and efficient neutron moderation. 




FIG. 1. DM-Icel7 detector schematic (left) and the location 
of the DM-Icel7 detectors within IceCube (right). A stainless 
steel pressure vessel encloses the Nal(Tl) crystal coupled to 
two PMTs, two IceCube mainboards, and two high voltage 
boards. The DM-Icel7 detectors are deployed 2457m deep 
in the South Pole ice and 7.5 m below the bottom IceCube 
modules. Det-1 is deployed below the DeepCore infill array 
(green volume), while Det-2 is located on the edge of the 1 km^ 
IceCube volume. 

DM-Icel7 has achieved an uptime of > 99% since the 
start of physics running in June 2011. This study uses 
3.5 yrs of data from July 1, 2011 to December 31, 2014. 
Day 1 is set to January 1, 2011; the data shown is from 
days 182 to 1460. 

Each DM-Icel7 detector contains a cylindrical 8.47 kg 
Nal(Tl) crystal (14.0 cm diameter, 15.0 cm height) cou¬ 
pled on each end to a photomultiplier tube (PMT). De¬ 
tails of the crystals’ history and contamination are pro¬ 
vided in [44] . The high voltage and digitizing electronics 
are located in an isolated volume above the upper PMT. 
The entire setup is enclosed in a stainless steel pressure 
vessel. 

Each DM-Icel7 PMT is independently controlled and 
monitored by its own IceCube mainboard (DOM-MB) 
assembly and high voltage (HV) board [50]. Waveforms 
are recorded and digitized when paired PMTs trigger 
within 350 ns of each other, referred to as local coinci¬ 
dence (“LC”). Timing is set by a clock on the surface and 
synchronized to the GPS receiver at the IceCube Lab, 
producing timestamps aligned with IceCube. Data is dig¬ 
itized into four channels to increase the dynamic range 
and read out across multiple time windows. The Flash 
Analog to Digital Converter (FADC) channel samples at 
40 MHz over 6.375 ps while the three 10-bit Analog Tran¬ 


sient Waveform Digitizer (ATWD) channels each record 
128 samples at a programmable rate (100 - 500 MHz) over 
a ^ 600 ns readout window that varies among the PMTs. 
The data presented here is from the ATWDO, ATWDl 
and ATWD2 channels, which have xl6, x2 and xO.25 
gain, respectively. The ATWD channels combine for a 
14-bit effective dynamic range, and events ranging from 
single photoelectrons to muons are recorded without sat¬ 
urating the data acquisition (DAQ). Single PMT thresh¬ 
old crossing rates are collected regularly as part of the 
monitoring data stream, which is recorded in addition to 
LC data collection. The trigger threshold for both LC 
and monitoring data is set to ^0.25 photoelectrons. De¬ 
tails of the DM-Icel7 hardware, deployment, and first 
data can be found in [44]. 

All figures shown are from Det-1. Det-2 exhibits con¬ 
sistent behavior and is included where it provides addi¬ 
tional information. 


III. MUON BACKGROUND 


Atmospheric muons are primarily produced via pion 
decay in the upper atmosphere. The muon production 
rate is dependent upon the likelihood of a pion to de¬ 
cay relative to an interaction with atmospheric nuclei. 
The mean free path of atmospheric pions depends on 
the atmospheric density, and changes in the density with 
temperature cause an annual modulation in the muon 
production rate with a maximum in the summer [42, 51]. 
Muons must have an energy of several hundred GeV in 
order to penetrate 2200m.w.e. of ice. Muon interactions 
are identified in DM-Icel7 by their characteristic pulse 
shape and large energy depositions (§111 A). The observed 
muon rate is 2.93 ± 0.04 p/crystal/day with a 12.3 ± 1.7% 
annual modulation (§111B). Muons interact as minimum- 
ionizing particles (MIPs) in the crystals (§111C). Muon 
identification is verified through coincidence with Ice- 
Cube (§111D), which is the subject of a separate pub¬ 
lication in preparation. 


A. Muon Event Identification 

Muon events in DM-Icel7 are separated from alpha 
and gamma events by a combination of pulse shape dis¬ 
crimination (PSD) and energy, as shown in Figure 2. 
Alpha interactions in Nal(Tl) produce pulses with a 
faster decay time than muon and gamma interactions 
[24, 52, 53], and they are identified through PSD, as 
shown in Figure 3. The PSD parameter, t, represents 
a charge-weighted mean time of the waveform. It is de¬ 
fined as: 
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where hi is the height of the waveform at time ti (de¬ 
scribed in detail in [44]). The r distributions are well 
separated, with less than one muon and one alpha event 
per year expected to be misidentified. 



Det-1 Pulse Height [a.u.] 


FIG. 2. Event classification in Det-1 from waveform 
heights and pulse shapes. Gamma events exhibit pulse 
heights < 325 (a.u.) and PSD r parameters > 177 ns; al¬ 
pha events comprise the lower PSD r band, reaching np to 
650 (a.u.) in pulse height. Muons make np the remaining 
events, with pulse heights > 325 (a.u.) and PSD r > 177 ns. 
The high energy pileup is due to DAQ channel satnration, 
with broad PSD values arising from PMT satnration. The 
high energy alpha population is comprised of pulses with mul¬ 
tiple events in a single readout (e.g., bismuth-polonium) that 
survive PSD removal. 



FIG. 3. Average waveforms, normalized in height, of muon 
(black), alpha (red), and gamma (blue) events in Det-1. Al¬ 
phas decay faster than gammas and muons, allowing their 
removal from the muon event sample. Lower energy muon 
events were chosen for this figure to avoid distortion from 
DAQ or PMT saturation. 

The largest energy depositions in the detector are from 
muon events (§ IIIC), and isolating events above the 
gamma spectral cutoff at 325 a.u. produces a highly pure 
muon event sample. This spectral cutoff is verified by 
GEANT4 simulation of the gamma backgrounds [44], and 
fitting the spectrum determines that fewer than 3 gamma 
events are expected in the muon sample every year. Be¬ 
cause the muon energy deposition is proportional to the 
path length through the crystal, muons that clip the edge 
of the detector deposit less energy and are indistinguish¬ 
able from gamma events. Muon path lengths through 
the crystal are simulated using the angular distribution 


of incident muons in order to determine the expected 
muon energy deposition. The shape of this expected dis¬ 
tribution indicates that <10% of muons appear below 
the gamma spectral cutoff and are misidentified as gam¬ 
mas. Reported muon rates correspond to those observed 
and do not attempt to correct for those muons that are 
misidentified as gammas. 

The ATWDl channel is used to identify muons be¬ 
cause of its higher gain. Muon energy analysis requires 
the ATWD2 channel due to ATWDl channel saturation 
above G.SMeVge (see §IIIC). Muon selection criteria use 
data from a single PMT for each detector because vari¬ 
ations in PMT response (e.g., gain, linearity, saturation 
energy) make single PMT analysis more effective in pro¬ 
viding separation. 


B. Muon Rate and Modulation 

The average muon flux is consistent between Det-1 and 
Det-2 and is measured to be 2.93 ± 0.04 p/crystal/day. 
The 2200m.w.e. ice overburden lowers the muon flux 
at DM-Icel7 by five orders of magnitude relative to the 
surface [54]. Muons lose roughly 400 GeV as they pass 
through the ice, and the muon flux above 1 GeV is sim¬ 
ulated to be 1.2 X 10“^ p/m^/s at this depth [55]. This 
simulated flux predicts the 14 cm diameter, 15 cm tall 
crystals to detect ~3 p/crystal/day, consistent with the 
observed rate. 

The muon flux through DM-Icel7, /^, is observed to 
modulate with a 12.3 ±1.7% modulation amplitude, as 
shown in Figure 4 and Table I. The rate is parameterized 
as an average muon flux, 7°, with a modulating compo¬ 
nent of amplitude A/^, period T, and phase to'. 

/^ = /°±AJ^,cos(^^(t-to)) (111.2) 

The modulation is expected to have a one year period 
and a maximum rate in the summer (i.e., January at the 
South Pole), although small deviations from a sinusoid 
are expected due to the imperfectly sinusoidal form of 
the temperature modulation and temperature variations 
from year to year. Fixing a one-year period yields a mod¬ 
ulation amplitude of A7^ = 0.31 ± 0.05 p/crystal/day 
with a phase of to = January 22 ± 9 days. When allowing 
the period to float, the best fit period is T = 371 ± 9 days. 

The modulation phase and period are consistent with 
those observed by IceCube [56], as shown in Figure 4. 
A quick rise to the maximum rate and slower decline 
are visible in the IceCube data; DM-Icel7 does not have 
enough statistics to observe this deviation from a si¬ 
nusoid. The modulation fits to a fixed period sinu¬ 
soid (x^/d.o.f. = 50.1/39) and a floating period sinu¬ 
soid (x^/d.o.f. = 49.6/38) are statistically indistinguish¬ 
able (< Icr separation). A scaled fit to the IceCube data 
with a floating relative modulation amplitude (y^/d.o.f. 
= 31.3/20) is preferred over a scaled fit with a fixed am- 
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plitude (x^/d.o.f. = 36.2/21) to greater than 2cr. All are 
preferred to greater than 5a over the null hypothesis. 



FIG. 4. Combined muon rate of Det-1 and Det-2 (black cir¬ 
cles, top figure, left axis), fit to a fixed-period sinusoid (best 
fit blue dashed line with la and 3a bands), and overlaid with 
IceCube rates passing an eight-fold simple majority trigger 
(SMT8) [56] (red circles, right axis). DM-Icel7 data is binned 
by month, while IceCube data is a daily average. Statistical 
error on the IceCube data shown is within the size of the 
data point. The modulation is consistent across both DM- 
Icel7 detectors and fits a 12.3 ± 1.7% modulation amplitude. 
It is consistent with both the sinusoidal fit and the IceCube 
data, and the residual of the fit to data is shown in the lower 
plot. 

The amplitude of the annual modulation in the atmo¬ 
spheric temperature and the corresponding muon rate 
are maximal at the South Pole due to its geographic lo¬ 
cation; daily temperature fluctuations do not mitigate 
the annual variation. The muon production mechanism 
induces a modulation amplitude that increases with en¬ 
ergy, leading to a depth dependent amplitude; higher 
energy muons increasingly match the fractional modu¬ 
lation of the temperature modulation, which is 11% at 
the South Pole [51]. The muon modulation has been 
measured by IceCube (1450-2450 m) at the South Pole 
to be 8.6±1.2% [42]. DM-Icel7’s location at the bot¬ 
tom of IceCube means that only higher energy muons 
reach it, leading to a larger modulation amplitude that 
corresponds more closely to that of the atmospheric tem¬ 
perature modulation [51]. 

C. Energy Deposition 

Simulations indicate that the muons reaching DM- 
Icel7 are most likely to exhibit a MIP energy loss as they 
pass through the crystal (~ 7 MeV/cm) with total energy 


deposition proportional to track length. Accounting for 
the predicted muon angular distribution and the crystal 
geometry, the muon path lengths range from 0 to 21cm 
with a strong peak at 15 cm. This corresponds to an ex¬ 
pected muon energy distribution peaked at 100 MeV with 
a flat continuum extending down to 0 MeV. 

The energy spectrum of Det-1 is shown in Figure 5 with 
events ascribed to gammas, alphas, and muons as de¬ 
tailed in §III A. A linear energy calibration is established 
from the locations of gamma peaks out to 2.6 MeVee [44]. 
Integrated charge in the ATWDl channel is used for the 
energy scale of alphas and gammas. At higher energies, 
the integrated charge no longer exhibits a monotonic rela¬ 
tionship with energy due to the effects of DAQ and PMT 
saturation; therefore, pulse height in the ATWD2 chan¬ 
nel is used for muons. Linear extrapolation of the gamma 
peak height calibration results in an observed muon spec¬ 
trum extending to only 23 MeVee, despite the expectation 
of a spectrum peaked at 100 MeV. Uncorrected PMT sat¬ 
uration effects are responsible for the difference. 
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FIG. 5. The energy spectrum of Det-1 gamma (blue), alpha 
(red), and muon events (black). Muons produce the highest 
energy events in the detector, extending out to 23 MeVee. De¬ 
tector response effects produce a reconstructed energy that is 
significantly smaller than the simulated energy deposition of 
~100MeV. The highest energy alpha peak (~ 7.5 MeVee) is 
from bismuth-polonium decays. 


D. IceCube Coincidence 

DM-Icel7 muon identification is verified through ob¬ 
served coincidence with IceCube. Det-1 is located in the 
center of the IceCube area, with full azimuthal accep¬ 
tance; Det-2 is located on the outer edge of IceCube, 
limiting its solid angle acceptance. Of muon events iden¬ 
tified in Det-1 over the 2012-2013 run year, 93% were 
coincident with muon events in IceCube; 43% of muon 
events in Det-2 were coincident over this same period. 
These coincidence rates are consistent with expectation 
once IceCube data selection efficiency is taken into ac- 
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TABLE I. The average rate (1°), modulation amplitude (A/^), phase (to), and best fit period (T) (see Eq. III.2) of the muons 
in DM-Icel7. The modulation amplitude is reported in [|r./crystal/day] and in [%]. The amplitude and phase are from a 
fixed one-year period modulation, with the best fit period reported in the right column. These values are consistent between 
detectors and agree with expectation. 



1° 

[p/crystal/day] 

Fixed Period Fit 

[p/crystal/day] [%] 

to 

[days] 

Floating Period Fit 

T 

[days] 

Det-1 

2.89±0.05 

0.39 ±0.07 

13.5 ±2.4% 

Jan 22 ±11 

387 ±16 

Det-2 

2.96±0.07 

0.33 ±0.07 

11.1 ±2.4% 

Jan 21 ± 12 

370 ± 14 

Combined 

2.93 ±0.04 

0.36 ±0.05 

12.3 ±1.7% 

Jan 22 ± 9 

371 ±9 


count. This agreement with expectation verifies the low 
leakage of gamma and alpha events into the muon sam¬ 
ple. 


IV. PHOSPHORESCENCE 

Nal(Tl) crystals emit photons primarily in fluorescent 
scintillation but also release photons from longer excita¬ 
tions, known collectively as phosphorescence or afterglow 
(§IV A). Phosphorescence is broadly classified as shorter- 
lived (<ms) and longer-lived (>ms) excitation. Muon 
events in DM-Icel7 induce long-lived phosphorescence, 
increasing in magnitude with muon energy. Comprised 
of up to tens millions of photons, these phosphorescent 
periods last for tens of seconds (§IVB). Deadtime in DM- 
Icel7 precludes a statement on short-lived phosphores¬ 
cence. 


A. Phosphorescence in Nal(Tl) 

Phosphorescence refers to Nal(Tl) scintillation longer 
than the ^250 ns primary fluorescent decay compo¬ 
nent [57]. A suite of Nal(Tl) experiments has observed 
phosphorescence with timescales ranging from 1.5 ps to 
45 days [24, 53, 58-68]. Differences in experimental re¬ 
sults indicate numerous phosphorescence decay channels. 
No correlations between phosphorescence decay time and 
crystal size or temperature have been observed in these 
studies. Phosphorescence has been observed following 
exposure to cosmic rays, radioactive sources, and ultra¬ 
violet (UV) and visible light. It is determined that the 
initial intensity of decays shorter than the second-scale is 
proportional to the intensity of incident radiation, while 
that of longer decays (minutes to hours) is proportional 
to the absorbed dose [58]. Both the ANAIS and KIMS 
dark matter experiments have measured elevated trig¬ 
ger rates following large energy depositions. ANAIS re¬ 
ports a 70 -100 ms [53, 59] and a minutes-long [68] decay, 
and KIMS reports a few seconds decay [24]. St. Gobain, 
the manufacturer which produced the DAMA crystals, 
has measured a phosphorescent component over hours to 
days following exposure to UV light [60]. By contrast. 


DAMA has not reported any observation of phosphores¬ 
cence in their analysis data sample. No details of a phos¬ 
phorescence study have been disclosed, but it has been 
reported that no events in the dark matter signal region 
are attributed to phosphorescence. Phosphorescence is 
generally comprised of low energy events, with the St. 
Gobain reporting phosphorescence reaching lOkeV [60]. 
For a detailed discussion of experimental results, see [69]. 

The precise mechanisms of longer (>ms) phosphores¬ 
cent decays are not entirely understood, but they are pos¬ 
tulated to include long-lived metastable activator states 
and traps from crystal defects and impurities, both in¬ 
trinsic and irradiation-induced. In the metastable inter¬ 
pretation, an electron excites the activator center to a 
state that is forbidden to decay to ground and requires 
additional energy, generally from thermal excitation, to 
reach a state that can decay. The decay time is then 
characteristic of the metastable state lifetime. In the trap 
interpretation, the charge carriers are trapped and ther¬ 
mally released to recombine; this interpretation can pre¬ 
dict trap depths and equilibrium conditions, as described 
in detail in [65]. Shorter phosphorescent decays are un¬ 
derstood as emerging from self-trapped exciton emission, 
lasting ps-ms, and a long-lived recombination channel 
that produces ps-long decays [70]. Long-lived recombi¬ 
nation occurs when an electron and hole reach different 
thallium (Tl) activator sites, and the electron must be 
thermally released before recombination [70, 71]. The 
phosphorescence intensity, spectral composition, and de¬ 
cay time are dependent upon the purity of the raw ma¬ 
terial, growing conditions, and irradiation exposure [70]. 
Phosphorescence is highly variable across crystals and 
must be studied individually for each detector. 


B. Muon-Induced Phosphorescence in DM-Icel7 

Muon interactions can induce phosphorescence in the 
DM-Icel7 crystals. Phosphorescence manifests as an in¬ 
crease in the number of uncorrelated single photons emit¬ 
ted directly following a muon event in DM-Icel7. In¬ 
creases of up to tens of millions of photons have been 
detected. Phosphorescence is observed in both LG data, 
which requires coincidence between the two PMT coupled 
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to the same crystal within 350 ns, and single PMT thresh¬ 
old crossings registered in monitoring data (see §11). 

The monitoring data stream, rather than the LC data 
stream, is used to determine phosphorescence decay pa¬ 
rameters because it is a more direct measurement of pho¬ 
ton emission rates. This analysis was performed on an 
eight month subset of the data (July 1, 2011-February 
28, 2012), during which monitoring data was recorded ev¬ 
ery 2s. Subsequent data was monitored every 60s [44], 
limiting its ability to study phosphorescence. The photon 
rate in the monitoring stream is determined by counting 
the number of threshold crossings over a 900 ms period. 
Photon rates for this study are determined by averaging 
the first two monitoring data records in each PMT after 
a muon event, which corresponds to averaging across the 
4 s following a muon interaction. 

A statistically signihcant event rate increase is induced 
by 96% of muons, which induce rates higher than ob¬ 
served in 90% of normal running distributions. The mag¬ 
nitude of phosphorescence increases with muon energy, as 
shown in Figure 6. Trigger rates during phosphorescence 
can spike outside of the fluctuations of average trigger 
rates, which are 2.5 Hz/crystal in coincident data and 
lOOHz/PMT in monitoring data. Single PMT threshold 
crossings in the monitoring stream have been observed to 
reach as high as 230 kHz during phosphorescence. The 
photon rate increases immediately following the muon 
events and decays exponentially with a 5.5 ± 0.5 s decay 
time, as shown in Figure 7. The brightest phosphores¬ 
cence events (<1% of phosphorescence) exhibit a second 
resolvable decay time of ~30s. 



FIG. 6. The average photon rate from the first two moni¬ 
toring data records following a muon event, averaged across 
both PMTs in Det-1 (black) and compared to the average 
photon rate (blue dashed), shown with single-sided bands en¬ 
compassing 90% (dark blue) and 99% (light blue) of data. 
The energy deposition is determined by the pulse height of 
the muon event in a.u. Data is from the eight month period 
when monitoring data was taken at 2 s intervals. Det-2 shows 
consistent behavior. 

The high rate of phosphorescent photon emission leads 
to accidental LC triggers during large phosphorescence 



FIG. 7. Muon-induced phosphorescence is observed as ele¬ 
vated photon rates that decay with a characteristic decay time 
of 5.5 ± 0.5 s. Each color corresponds to a separate phospho¬ 
rescence period in Det-1 monitoring data. The x-axis is set to 
zero at the time of the muon event inducing phosphorescence. 
The largest phosphorescence periods, including the black and 
red curves shown, exhibit a second decay time of ~30 s. Det-2 
shows consistent behavior. 

periods. Phosphorescence induces statistically significant 
LC rate increases following 2% of muons, when they are 
identified in LC data by recording the number of trig¬ 
gers in the 30 s after the muon and comparing this to 
a baseline rate defined as the number of triggers in the 
30 s before the muon. The elevated LC trigger rates dur¬ 
ing phosphorescence exceed the designed range for the 
DAQ. The 730 ps deadtime that follows each event [72] 
and the scale of observed phosphorescence prevents the 
detection of the ms-scale decays reported by other ex¬ 
periments [53, 59]. An additional deadtime is introduced 
above a LC rate of 50 Hz, which only occurs during phos¬ 
phorescence. Analysis of these considerations and inclu¬ 
sion of quantum efficiency indicates that only a fraction 
of photons are recorded in most phosphorescence events, 
falling with increasing trigger rate to only 0.03% of pho¬ 
tons recorded in the largest, DAQ saturated events. Of 
the ~10® photons emitted in the largest phosphorescent 
event, 3940 triggers were recorded. When corrections 
for these effects are implemented, the photon rate ex¬ 
trapolated from LC data is consistent with the recorded 
monitoring rate. LC phosphorescence data from the en¬ 
tire 3.5 yr dataset is used for analysis, and it is consistent 
before and after the monitoring frequency change. 

Phosphorescence events are primarily comprised of dis¬ 
crete, low energy pulses that are dominated by sin¬ 
gle photoelectron-like triggers and are rejected as noise. 
Det-1 (Det-2) has a light yield of 5.9 (4.3) photoelec- 
trons/keV, and energy is defined as the integrated charge 
over the ~600 ns readout window. Figure 8 shows the 
energy spectrum of all LC phosphorescence triggers in 
Det-1, isolated by subtracting a baseline spectrum from 
the spectrum of all pulses in the 30 s after the muon. 
The baseline spectrum is made from all pulses in the 30 s 
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before the muon. LC phosphorescence triggers produce 
a spectrum that is peaked below 1 keVee with no statis¬ 
tically significant component above 2 keVee even before 
noise rejection. Phosphorescence events do not mimic 
normal scintillation events and are therefore rejected by 
the standard noise removal procedures described in [44]. 
This rejection leaves a spectrum that is consistent with 
the average, non-phosphorescence background spectrum 
at all energies. Muon-induced phosphorescence therefore 
does not impact the dark matter analysis sample in DM- 
Icel7. 



FIG. 8. ATWDO energy spectrum of phosphorescent LC trig¬ 
gers in Det-1 (black), isolated by subtracting the baseline 
spectrum derived from all pulses in the 30 s before the muon 
(blue) from the spectrum of the 30s after the muon (red). 
Muons that saturate the DAQ or are within 30 s of a run 
transition have not been included. The residual indicates the 
additional events due to phosphorescence, which appear as 
low energy noise. The inset shows the spectral comparison, 
with statistical error bars, after the application of noise cuts 
and illustrates the rejection of phosphorescence events with 
noise removal. 

Phosphorescence is consistent across both PMTs on 
each crystal and across PMT high voltage settings in DM- 
Icel7. This indicates that it is not the result of an effect 
from the PMT dynode chain. 

V. IMPLICATIONS FOR SODIUM IODIDE 
EXPERIMENTS 

Phosphorescence induced by muons has been proposed 
as a potential modulating background in the region of 
interest for Nal(Tl) dark matter experiments [29-31]. 
Such an effect requires a sufficient photon rate follow¬ 
ing a muon to produce photoelectron pileup that re¬ 
constructs above 2keVee- Muon-induced phosphores¬ 
cence is not rejected by traditional muon tagging and 
can evade multi-crystal and energy cuts. This must be 
considered when assessing the feasibility of muons and 
muon-induced phosphorescence to mimic a low energy 
scintillation-like signal. 


The observed phosphorescence in DM-Icel7 indi¬ 
cates that this is unlikely to be the sole source of 
the annual modulation signal observed in DAMA/Nal 
and DAMA/LIBRA, consistent with those experiments’ 
statements [32-34], although variations in detector de¬ 
tails and in phosphorescence across crystals complicate 
a comparison between experiments. In order for phos¬ 
phorescence to be the sole source of DAMA signals, the 
events in the region of interest must be entirely muon- 
induced, and this is highly improbable [32, 33]. Applying 
the phosphorescence observed in DM-Icel7 to DAMA’s 
muon background indicates that the event rate would 
be too small to produce the observed modulation, even 
before noise removal. Production of the DAMA modula¬ 
tion signal would require a significantly higher number 
of phosphorescent photons with a time structure con¬ 
ducive to photoelectron pileup. Additionally, the phase 
of the modulation would have to precede that of muons 
by 45 days [9, 43]. A second annually modulating back¬ 
ground with a different phase could alter the apparent 
amplitude and phase of the combined signal [30], but no 
such background has been established [34, 36]. 

DAMA has reported that muons and associated phos¬ 
phorescence cannot produce their observed annual mod¬ 
ulation signal [32-34]. The observations of DM-Icel7 
presented here provide no indication that muon-induced 
phosphorescence is responsible for the modulating signal 
observed in DAMA/Nal and DAMA/LIBRA. However, 
variations in the properties of phosphorescence among 
crystals [59, 70] require that each Nal(Tl) experiment 
study the possible signal leakage of its phosphorescence 
background. 


VI. CONCLUSIONS 

We have identified and characterized muons and 
muon-induced phosphorescence in DM-IceI7. Muon 
events are identified by their pulse shape and en¬ 
ergy deposition, yielding an average observed flux of 
2.93 ±0.04 p/crystal/day. The muon rate, correlated to 
atmospheric temperature, modulates with a 12.3 ±1.7% 
amplitude and a phase corresponding to a maximum flux 
on January 22 ±9 days. Muon identification is validated 
with observed coincidence in IceCube, indicating a highly 
pure muon sample in DM-Icel7. Muon events induce a 
long-lived phosphorescence that is comprised of up to 
tens of millions of photons. The rate of photons from 
phosphorescence falls off exponentially with a 5.5 ± 0.5 s 
decay time. While this effect provides a mechanism for 
muons to induce a modulation at low energies that es¬ 
capes muon removal cuts, the phosphorescence pulses 
observed in DM-Icel7 appear primarily below 2keVee- 
These events are rejected as noise, and they do not have 
the event rate in DM-Icel7 required to produce a mod¬ 
ulation consistent with the DAMA signal, even before 
noise removal- Variations in detector details and in the 
properties of phosphorescence between crystals render an 












inter-experimental comparison difficult and motivate de¬ 
tailed studies by all experiments on the potential impact 
of phosphorescence in the region of interest. We con¬ 
clude, however, that the annually-modulating signal re¬ 
ported by the DAMA experiments cannot be accounted 
for by phosphorescence with the specihc characteristics 
observed in DM-Icel7. 
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